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BENDING WAVES AND ORBITAL INCLINATIONS. 

Wm. R. Ward, Jet Propulsion Laboratory, California Institute of Technology, 
Pasadena, CA 91109, and San Juan Capistrano Research Institute, San Juan 
Capistrano , CA 92675 . 

Disc tides may play an important role in the formation of a planetary 
system. Modification of protoplanet semi -major axes and eccentricities 
through density waves was first suggested by Goldreich and Tremaine (1979, 
1980) and calculations of linear and non-linear protoplanet-nebula 
interactions have been carried out by a number of researchers (e.g., 
Goldreich and Tremaine, 1979; Papaloizou and Lin, 1984; Lin, Papaloizou, and 
Savonije , 1990; and Ward, 1986, 1990). The possible significance of radial 
drift to the accretion process has been discussed by Hourigan and Ward 
(1984) and Ward (1986, 1989a), while the evolution of orbital eccentricities 
has been studied by Goldreich and Tremaine (1980, 1981) and Ward (1989b) for 
a perturber orbiting external to a ring and by Ward (1988) for a perturber 
embedded in the ring. Here we stress the possible importance of bending 
waves to the early evolution of protoplanet inclinations. 

Resonant interactions between a disc and a perturber in an inclined 
orbit occur at locations of vertical resonance, which are analogues of 
Lindblad resonances and for which the disc disturbance takes the form of a 
spiral bending wave (e.g., Shu et al., 1983; Shu, 1984). The perturbing 
vertical acceleration can be Fourier decomposed into terms of the form g p - 
Re{/ m exp i(m0-wt)}. The strongest terms have frequencies - (ra-e)O p and 
amplitudes 
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where r s , O g , and M # are the semi-major axis, mean motion, and mass of the 
perturber, I is its inclination, and b“ / 2 (7) is a Laplace coefficient with 
argument 7 - r/r # . For each m there are two forcing terms given by e - ±1. 
Resonance occurs when the vertical oscillation frequency, O z , equals the 
Doppler shifted forcing frequency, i.e., when 
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where 0 is the local disc orbital frequency. In a Keplerian disc about a 
primary, M p ; Q - 0 ^ - (GM p /r 3 ) 1/2 , and resonances fall at 7 ~ l+4e/3m. 

Borderies et al. (1984) determined the effects of the reaction torque 
from an m th order resonance on the inclination of the perturber and found 
that it increased at a rate 



where d ■ rdD/dr, and a is the surface density of the disc. Both inner and 
outer resonances contribute to excitation. For a Keplerian disc, l®l - 
Slm+eln 2 . If |1 -t| « 1, Laplace coefficients can be approximated by 
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modified Bessel functions, b£ /2 ( 7 ) * (3m 2 ^w)^ (m| 1- 7 | ) . These are to be 
evaluated at resonance, i.e., m|l- 7y | - 4/3. For m » 1, eqn (3) becomes, 
(esc I)dl/dt - 0.0235(M s /M )0(<rr 2 /M p )m 3 . To find the total change in I, eqn 
(3) must be summed over all resonances. 

As an application, consider a Jovian size protoplanet forming from a 
circumstellar disc. The protoplanet is assumed to have terminated its gas 
accretion phase by tidally truncating the disc as proposed by Lin and 
Papaloizou (1979) . Truncation will abort radial drift and damp the 
eccentricity, provided corotation torques from the edges of the gap dominate 
Lindblad torques (e.g., Goldreich and Tremaine, 1980, 1981; Ward, 1989b). 
However, Borderies et al. (1984) point out that there are no analogous 
corotation torques in the vertical resonance problem and thus no analogous 
mechanism to halt inclination growth. Resonances are summed to m — 
(4/3) (r/w), where w is the gap half-width. The resulting rate implies a 
characteristic inclination growth time, r - sin I/(dI/dt) , of 
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The minimum size gap will be on the order of the scale height of the disc, h 

- c /fi ~ .025r(aT 2 ) 1/2 , where a is in AU and T 2 is the disc temperature in 

10 2 °K. Note that perturbation wavelengths, A - 2nr/m £ h, so that the thin 
disc model used to derive eqn (3) should be roughly applicable. Adopting a 
one solar mass primary, M p - 2 X 10 33 gm. , a Jovian size protoplanet, M s /M 
-10 , at 5.2 A.U. , and a nebula surface density and temperature of , *4o6 

gm/cm 2 , and 160°K, respectively; eqn (4) reads, r - 10 3 (w/h) A years. The 
expected lifetime of the nebula is of order 10 7 years (Adams and Shu, 1986; 
Walter, 1986). A characteristic growth time this long would require a gap 
width of order 10 scale heights. 

The critical protoplanetary mass necessary to open a gap against 
viscous diffusion depends on the turbulent viscosity, u , of the disc, i.e., 
M c /M p - (iVr 2 n) 1/2 (c/rfl) 3/2 (e.g., Hourigan and Ward, 1984; Papaloizou and 
Lin, 1984). The mass required to sustain the edge at a distance w > h, 
exceeds M e by a factor (w/h) 3/2 . Thus, a gap width w - lOh could be 
maintained by a protoplanet £ 30 times that needed to truncate the disc. 
However, since disc truncation will inhibit further gas accretion, it is not 
obvious how to reconcile these requirements. Borderies et al. (1984) 
speculated that the sign of dl/dt may change when I t w/a, citing a similar 
behavior in the eccentricity case found by Goldreich and Tremaine (1981). 
Nevertheless, a minimum width gap would still permit inclination growth to I 

- h/r « 4° for both Jupiter and Saturn, which seems inconsistent with their 
small mutual inclination. Alternatively, eqn (4) could imply that much of 
the gas accretion phase of the giant planets was correlated with disc 
dispersal so that their co-existence was for an interval significantly 
shorter than the disc's lifetime. 

We should point out that the excitation rate of eqn (3) was calculated 
for a thin, self- gravitating particle disc. By contrast, a minimum mass 

solar nebula is a pressure dominated disc with a finite thickness in 
hydrostatic equilibrium with the vertical component of the solar gravity. 
However, in the case of density waves, it is well known that the same 
angular momentum flux is carried by both gravity and pressure waves and that 
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the torque at a Lindblad resonance is independent of the restoring force 
responsible for the wave action. We suspect this is true for bending waves 
as well and that the torque strength will prove rather insensitive to the 
details of the wave process. The application discussed here underscores the 
importance of developing a more rigorous treatment of this problem. 

REFERENCES 

Adams, F. C., and Shu, F. H. (1986). Astrophys . J. , 308, 836. 

Borderies, N. P., Goldreich, P., and Tremaine, S. (1984). Astrophys . J. , 

284, 429. 

Goldreich, P, and Tremaine, S. (1979). Astrophys. J., 233, 857. 

. (1980). Astrophys. J., 241, 425. 

. (1981). Astrophys. J., 243, 1062. 

Hourigan, K. and Ward, W. R. (1984). Icarus, 60, 29. 

Lin, D. N. C., and Papaloizou, J. (1979). M.N.R.A.S., 186, 789. 

Lin, D. N. C. , Papaloizou, J., and Savonije, G. (1990). Astrophys. J. , 364, 
326. 

Papaloizou, J. and Lin, D. N. C. (1984). Astrophys . J., 285, 818. 

Shu, F. H. , Cuzzi, J. N. , and Lissauer, J. J. (1983). Icarus, 53, 185. 

Shu, F. H. (1984). "Waves in planetary rings", in Planetary Rings, 

(Greenberg and Brahic, Eds.), University of Arizona Press, pp. 513-561. 
Walter, F. M. (1986). Astrophys. J., 306, 573. 

Ward, W. R. (1986). Icarus, 67, 164. 

. (1988). Icarus, 73, 330. 

. (1989a). Astrophys. J., 345, L99. 

. (1989b). Astrophys. J., 326, 526. 

. (1990). in preparation. 


